A comment on AdS collapse of a scalar field in higher dimensions 
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We point out that the weakly turbulent instability of anti-de Sitter space, recently found in [3| 
for four dimensional spherically symmetric Einstein-massless-scalar field equations with negative 
cosmological constant, is present in all dimensions d + 1 for d > 3, contrary to a claim made in Q]. 
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In a recent paper [H two of us reported on numeri- 
cal simulations which indicate that anti-de Sitter (AdS) 
space is unstable against the formation of a black hole 
under arbitrarily small generic perturbations. This insta- 
bility was conjectured to be triggered by a resonant mode 
mixing which moves energy from low to high frequencies. 
Although the simulations of [1] were done only in four 
spacetime dimensions, it was clear from the nonlinear 
perturbation analysis given there that the same weakly 
turbulent mechanism of generating instability operates 
in all supercritical dimensions d -f- 1 for d > 3. Thus, 
we found it surprising that Garfinkle and Pando Zayas, 
who later looked at the same problem in 4 + 1 dimen- 
sions (apparently unaware of Ref.fr] as they did not cite 
it), wrote in [2]: "We [...] establish that for small val- 
ues of the initial amplitude of the scalar field there is no 
[sic!] black hole formation, rather, the scalar field per- 
forms an oscillatory motion typical of geodesies in AdS." 
The purpose of this comment is to show that, contrary 
to the above quoted claim, AdSd+i is unstable against 
gravitational collapse for all d > 3. To this end, we first 
generalize the formalism of T] to higher dimensions and 
recall the key argument for collapse of arbitrarily small 
generic initial data. Second, using our code in d = 4 we 
evolve numerically the same initial data as those of 0| 
and show that, as expected, they collapse. Finally, we 
try to identify a possible source of the error in Q . 

We parametrize the (d-f l)-dimensional asymptotically 
AdS metric by the ansatz 

ds^ = — — (~Ae-^^dt^ + A-^dx^ + sin^x dnj,) , 
cos^x ^ ' 

(1) 

where = — d(d — 1)/(2A), dn^_^ is the round metric 
on —oo <t<oo, 0<x< 7r/2, and A, 6 are 

functions of (t, x). For this ansatz the evolution of a self- 
gravitating massless scalar field (f>{t, x) is governed by the 
following system (using units in which SttG = d — 1) 

$=(Ae"^n)', rr=_L_(tan'^-i.Ae-^$)', (2) 



d - 2 + 2sin^x 

1' = (1 

sm x cos X 

5' = - sin a: cos X -I- n^) , 



— - ' - ~ A)- sinzcosxA ($2 +n2) , 

sin X cos X 

(3) 

s2 , tt2\ ^^-j 



where ' = dt, ' — dx, — <f)' and 11 — A^^e^cj). We want 
to solve the system (2-4) for small smooth perturbations 
of AdS solution 4> = A — 1,5 = Q. Smoothness at the 
center implies that near x = Q 



^{t,x) = fo{t) + 0{x'), 
A{t,x) = l + 0{x^) , 



Sit,x) = Oix^), 



(5) 



where we used normalization 6(t,0) = so that t is the 
proper time at the center. Smoothness at spatial infinity 
and finiteness of the total mass M imply that near x = 
tt/2 we must have (using p = tt /2 — x) 



Cf>{t,x)^foo{t)p'' + O{p'+^) 

A{t, x) = l~ Mp'^ + O {p'^+^) 



5{t,x)^5^{i) + 0{p'^) 
(6) 



where the power series expansions are uniquely deter- 
mined by M and the functions /oo(i), 5oo{t) (which in 
turn are determined by the evolution of initial data). One 
can show that the initial-boundary value problem for the 
system (2-4) together with the regularity conditions (5) 
and (6) is locally well-posed. 

In []| the instability of AdS^ was conjectured to result 
from the resonant mode mixing which moves energy from 
low to high frequencies. It was argued that this process of 
energy concentration on increasingly small spatial scales 
must be eventually cut off by the formation of a black 
hole. It is easy to see that the same mechanism is at 
work for all d > 3. This follows from the fact that, using 
the PDE terminology, the system (2-4) is fully resonant. 
More precisely, the spectrum of the linear self-adjoint 
operator, which governs the evolution of linearized per- 
turbations of AdSd+^, L ~ — tan^"''' x dx {tan^~^ x dx) , 

2j)2, {j = 0,1,...). The key 



is given by 



(d 



point is that the frequencies ojj are equally spaced, so 
already at the third order of nonlinear perturbation anal- 
ysis one gets resonant terms for any frequency ujj such 
that j = ji + j2 — js , where jk are indices of eigenmodes 
present in the initial data. Some of these resonances lead 
to secular terms which signal the onset of instability at 
time t = 0{e^'^), where e measures the size of initial 
data. 

To substantiate this heuristic argument we solve the 
system in d = 4 numerically, using the method of 

[l|. For easy comparison of results, we take the same 
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approximately ingoing Gaussian initial data as Ref.Q]: 
$(0, x) = 9^,0(0, x) = n(0, x), where 

^ e / (tanx-ro)2\ 
0(0, a;) = exp 1^ — j (7) 

with — A, a — 1.5. As in Q we set the AdS radius 
£ = 1, hence their and our radial coordinates are related 
by r = tanx, while the time coordinates are identical. 
We denote their amplitude A by e (to avoid conflict of 
notation); the factor comes from the difference in 

units: we use SttG = 3, while in Q SttG = 1. Note 
that the data Q slightly violate the regularity condition 
([5]) since $(0,0) is not exactly zero, however an error 
generated by this 'corner singularity' is negligible. 
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FIG. 1; Time of horizon formation tn vs amplitude for ini- 
tial data (O. The upper plot depicts the first 'step' of the 
'staircase' function tH{e) corresponding to large data solu- 
tions which collapse during the first implosion. The horizon 
radius varies from xh = 0.5 to zero (from right to left). For 
£ > 0.005, the plot coincides with Fig. 4 of 01 which, veri- 
fies that our results agree with those of for short enough 
times. The lower plot depicts a few further 'steps' of tnis) 
corresponding to solutions which bounce several times from 
the AdS boundary before collapsing. 

For large e, the solution quickly collapses (the forma- 
tion of an apparent horizon is detected by the metric 
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FIG. 2: Upper plot: the upper envelope of n^(f, 0) for initial 
data (O with three relatively small amplitudes. After mak- 
ing 95 (for e = 0.0002), 198 (for e = ^2 • 0.0001, and 405 
(for e = 0.0001) reflections, all solutions eventually collapse. 
Lower plot: the curves from the upper plot after rescaling 
£~^n^(£^t, 0) seem to converge to a limiting curve. 

function A{t,x) touching zero at some xh)- As e de- 
creases, the horizon radius takes the form of the right 
continuous sawtooth curve xh{£) (see Fig. 1 of [1]) with 
jumps at critical points £„ where lim^^^+ xnis) — (the 
index n counts the number of reflections from the AdS 
boundary before collapse) . Accordingly, the time of hori- 
zon formation i/f(e) is a monotone decreasing piecewise 
continuous function with jumps at each e„ (see Fig. 1). 

For small initial data the weakly nonlinear perturba- 
tion analysis described in [1| predicts the onset of insta- 
bility at time t ^ e^^. Numerics indicates that for suf- 
ficiently small e this scaling holds approximately almost 
all the way to the collapse, that is i//(e) ~ e~^. The 
evidence for this fact is shown in Fig. 2 which depicts the 
evolution of three solutions with small amplitudes differ- 
ing by a factor of v^. Note that this scaling implies that 
the computational cost of numerical evolution increases 
rapidly as e decreases (since solutions have to be evolved 
for longer and longer times on finer and finer grids). 
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Finally, let us make a few remarks about the paper . 
The content of that paper, when stripped of the standard 
promotional material for the application of AdS/CFT 
correspondence to the quark-gluon plasma, boils down 
to two statements concerning the solutions of five di- 
mensional spherically symmetric Einstein-massless-scalar 
field equations with negative cosmological constant: (i) 
large initial data lead to collapse and the larger the data, 
the shorter the time of horizon formation; (ii) small ini- 
tial data do not form black holes. The statement (i) is 
true but trivial. The statement (ii), as we have shown 
above, is false. We wonder what led the authors of 0] to 
reach this conclusion. The numerical method used in Q 
is rather crude; it is based on the second order finite dif- 
ference code on a fixed nonuniform grid. The radial coor- 
dinate r is not compactifled and the AdS timelike bound- 
ary at r = oo is mimicked by an artificial reflecting mirror 
at r,nax = 10. Although this fact might make long-time 
quantitative results somewhat inaccurate (since the pulse 
gets reflected before reaching the true boundary), we do 
not think it had any effect on the claim (ii). We think 
that the culprit was the lack of sufficient resolution. In 
the first (vl) version of [2|, Fig. 2vl shows three implo- 
sions for the solution with amplitude e = 0.001, which 
(together with the discussion in the text) seems to sug- 
gest implicitly that no black hole forms later and "the 
scalar field performs an oscillatory motion" forever. We 
evolved these data with our code and observed horizon 
formation at xh — 1-88 • 10~^ during the first implo- 
sion at time tn = 1.398 which means that the absence of 
collapse and the "oscillatory motion" of Fig. 2vl are nu- 
merical artifacts. It appears that the authors of [2| have 
realized that their numerical simulations suffered from 
insufficient spatial resolution because in the second (v2) 
version of |2| the number of grid points was increased 
from 800 to 6400. We gather that this improvement of 
resolution has helped to detect the collapse in Fig.2vl 



during the first implosion because the new Fig. 2v2 de- 
picts four implosions of a 'non-collapsing' solution with 
five times smaller amplitude e = 0.0002. But this has 
just postponed the problem of the loss of resolution until 
a later time (not shown in the plot). Indeed, our numer- 
ical evolution of these data (shown in Fig. 2) yields hori- 
zon formation at xh ~ 4.8 • 10~^ after time tn ~ 299.8. 
Even if the numerical simulation used to produce Fig.2v2 
had been run for such a long time, it would not have had 
enough resolution to capture the spatio-temporal struc- 
ture of horizon formation since this structure develops 
below the first point of their grid. For comparison, our 
fourth-order code with global adaptive mesh refinement 
reached the level of 2^ ''-I- 1 grid points just before collapse. 
Anyway, the lesson is that the long-time numerical simu- 
lations of asymptotically AdS spacetimes are challenging 
even in spherical symmetry and one should be careful 
in jumping to conclusions about the late time dynam- 
ics, especially without an analytic understanding of the 
problem. 
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